UNIQUENESS FOR INVERSE STURM LIOUVILLE PROBLEMS 
WITH A FINITE NUMBER OF TRANSMISSION CONDITIONS 
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Abstract. We establish various uniqueness results for inverse spectral prob- 
lems of Sturm— Liouville operators with a finite number of discontinuities at 
interior points at which we impose the usual transmission conditions. We 
consider both the cases of classical Robin and of eigenparameter dependent 
boundary conditions. 
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1. Introduction 

In the seminal paper [7], Hald, motivated by the inverse problem for the torsional 
modes of the earth, investigated Sturm-Liouville problems with a discontinuity at 
an interior point. Hald proved a Hochstadt-Liebermann result in the case of one 
transmission condition which was later on extended to two transmission condi- 
tions by Willis [19]. Moreover, Kobayashi [9] proved a similar result in the case 
for problems with a reflection symmetry. More recently, Mukhtarov, Kadakal and 
Muhtarov [12] and two of us [15] have investigated the case with one transmission 
condition and eigenparameter dependent boundary conditions, and derived asymp- 
totic formulas for the eigenvalues and eigenfunctions. Even more recently, these 
results were extended to two and three transmission conditions in [8] and [14) . 
respectively. The purpose of the present paper is to show how to handle an arbi- 
trary finite number of transmission conditions and to use the asymptotic formulas 
to prove several uniqueness results. In particular, we will introduce a Weyl m- 
function which uniquely determines the parameters of the problem. We also show 
that this Weyl function is a meromorphic Herglotz-Nevanlinna function which is 
uniquely determined by its poles and residues, as well as by its poles and zeros. 
In particular, we also obtain a two spectra result. This generalizes the results of 
Amirov [3] in the case of one transmission condition to the case of a finite number 
of transmission and eigenparameter dependent boundary conditions. Moreover, we 
will also generalize the Hochstadt-Liebermann type result from Hald to the present 
situation. 

To the best of our knowledge, this is the first result concerning more than three 
transmission conditions. In particular, it was necessary to modify the usual argu- 
ments at several places in order to make up for some key estimates which cannot 
be easily shown in the present situation (cf. the intricate nature of the high energy 
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asymptotics of solutions in Theorem 13. 1|) . For related results, we refer to [I], [5], 

m, usi, m, p. 

Sturm-Liouville problems with transmission conditions at interior points arise 
in a variety of applications in engineering and we refers to 3 for a nice discussion 
and further information. Here we only want to mention that they also appear 
in the description of delta interactions (which play an important role in quantum 
mechanics 2.) and of radially symmetric quantum trees (cf. the discussion in 
Section 4 of [T3] and the references therein). For general background on inverse 
Sturm-Liouville problems we refer (e.g.) to the monographs [5J, |llj . |17j . 

We will first start with the usual Robin boundary conditions and then briefly 
show how to extend the present approach to the more general case of eigenparameter 
dependent boundary conditions in our last section. 

2. The Hilbert space formulation and properties of the spectrum 
In the first part of our paper we consider the boundary value problem 

(2.1) £y := -y" + qy = Xy 
subject to the Robin boundary conditions 

Li{y):=y'(0) + hy(0) = 0, 

(2.2) L 2 (y):=y'(ir)+Hy(ir)=0 

with transmission (discontinuous) conditions 

Ui{y) := y(d t + 0) - cny(di - 0) = 0, 

(2.3) V t (y) := y'{d, + 0) - &#'(* - 0) - ay fa - 0) = 0, 

where q(x) is real- valued function in L [0, 7r]. We also assume that h, H and a%, 
bi, Ci di, i — 1,2, . ..,m — 1 (with m > 2) are real numbers, satisfying aibi > 0, 
da = < d\ < d 2 < ... < d m — i < d m = n. For simplicity we use the notation 
L = L(q(x); h; H; d t ), for the problem (J2~T]) ~(|23 | . 

To obtain a self-adjoint operator we introduce the following weight function 

'1, 0<a:<di, 



(2.4) w(x) 



i 

aibi ' 



d\ < x < d 2 , 



< x < IT. 



. a\b\ ■■■a m _ib TT1 _i ' 

Now our Hilbert space will be T-L := ^2((0, 7r); w) associated with the weighted inner 
product 



(2-5) (f,g) n := / fgw. 

Jo 

1/2 

The corresponding norm will be denoted by ||/||« = (/, f)q[ ■ In this Hilbert space 
we construct the operator 

(2.6) A-.n^n 

with domain 

fJ'eACiu™- 1 ^,^^)), 



(2.7) dom(A) = ifen 



^eL 2 (0,7r), Ui(f) =Vi{f)=0 
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by 

Af = if with/ e dom {A) . 

Throughout this paper AC{ U™ _1 (di, dj+i)) denotes the set of all functions whose 
restriction to (di, d»+i) is absolutely continuous for all i = 0, . . . , m — 1. In partic- 
ular, those functions will have limits at the boundary points di. 

Lemma 2.1. The operator A is self-adjoint. 

In particular, the eigenvalues of A, and hence of L, are real and simple. To 
see that they are simple it suffices to observe that the associated Cauchy problem 
([2"Tjl . (|2~3|) subject to the initial conditions f(x ± 0) = f , f(x ± 0) = f x (with 
.To G [0, 1]) has a unique solution. 

For any function / £ dom (A) we will denote by fj, 1 < j < to, the restriction of 
/ to the subinterval (dj_i,dj). Moreover, we will set fj{dj-i) = f(dj-i +0) and 
fj(dj) = f(d j -0). 

Suppose that the functions </?(a;, A) and ip(x, A) are solutions of (|2.1|) under the 
initial conditions 

(2.8) ^(0,A) = 1, ^'(0,A) = -/ l; 
and 

(2.9) ^(7r,A) = l, ip'(n,\) = -H 

as well as the jump conditions p.3[) . respectively. It is easy to see that equation 
(|2.ip under the initial conditions (|2.8j) or (|2.9|) has a unique solution yi(a;, A) or 
ip m (x 1 X), which is an entire function of A G C for each fixed point x e [0,di) or 
x G (d m _i,7r]. From the linear differential equations we obtain that the modified 
Wronskian 

(2.10) W(u,v) = w(x)(u(x)v'(x) - u'(x)v(x)) 

is constant onie [0, di) U™ -2 (di,di + 1) U (d m _i,7r] for two solutions £u = Aw, 
£v = Xv satisfying the transmission conditions (|2.3p . Moreover, we set 

(2.11) A(A) := T4%(A), V(A)) = £i(V(A)) = -w(7r)L 2 (^(A)). 

Then A(A) is an entire function whose roots A„ coincide with the eigenvalues of 
L. Moreover, the eigenfunctions pi(x, X n ) and ipi(x,X n ) associated with a certain 
eigenvalue A n , satisfy the relation ipi(x, A„) = /3„ipi(x, A„), where, by 



(2.12) /9n = '0(O,A„). 

We also define the norming constant by 

7„ := \\(f(x, X n )\\ n 2 . 
Then it is straightforward to verify: 
Lemma 2.2. All zeros X n of A(A) are simple and the derivative is given by 

(2.13) A(A„) = -i^Pn- 

Finally, we point out a simple unitary transformation for our eigenvalue problem 
which is easy to check: 
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Lemma 2.3. The map 



maps A unitarily to A associated with &i = (ai/bi) 1 ' 2 , b\ = (bi/ai) 1 ' 2 , c\ — 
Ci(dibi)~ 1 ' 2 and all remaining items unchanged. In particular, dibi = 1 and hence 
w(x) = 1. 

Remark 2.4. After a similar transformation as above we can assume that a% = 1 
without loss of generality and then our operator is a special case of a measure-valued 
Sturm- Liouville operator [4j 



^ = W)i(- w{xW{x)+ Io v{t)dx{t) ) 



associated with the measure-valued potential 

m — l 

dx(x) = q(x)dx + J^ w(di+)ciSdi(x), 



where 6j is the Dirac delta 



located at d. 



3. Asymptotic form of solutions and eigenvalues 

Theorem 3.1. Let X = p 2 and t :— Imp. For equation (|2.1j) with boundary 
conditions (|2.2[) and jump conditions (|2.3[) as |A| — > oo, the following asymptotic 
formulas hold: 



(3.1) 



ip(x,X) = < 



cos px + O C* m p Tlx > ), 0<x<d x , 

ct\ cospx + a[ cosp(x — 2d\) + 0( cxp ^' T ' x ' ), di < x < d 2 , 

ct\u 2 cos px + cJ x a 2 cos p(x — 2d\) + a\cJ 2 cosp(x — 2d 2 ) 

+a[a 2 cosp(x + 2di - 2d 2 ) + 0{ cxpiMx) ), d 2 < x < d 3 , 



aia 2 ...a m -i cospx+ 

+a' 1 a 2 ...a m -i cos p(x — 2d{) + ... 
+aia 2 ...a / m ^ 1 cos p(x - 2d m _i)+ 
+a[a 2 a 3 ...a m -i cos p(x + 2d x - 2d 2 ) + .. 
+ai...a' i ...a'y..a m ^i cos p(x + 2di — 2dj) 



-ai...a' i ...a'y..a' k ...a m -i cos p(x — 2di + 2dj — 2dj~) 
-ai4...a^_! cosp(x + 2(~l) m - 1 d 1 + 2( 



-l) m - 2 d 2 



2d 



m-l) 



-o{ 



cxp(|-r|3j) 



), 



d m -l < X < 7T, 
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and 

(3.2) 



ip'(x,X) 



' p[— sinpx] + 0(exp(|T|x)), < x < d\, 

p[— a\ svapx — a[ sinp(:r — 2di)] + 0(exp(|r|x)), d\ < x < cfe, 

p[— a\a 2 sin pa; — a' x a% sinp(x — 2d\) — 

— a\a' 2 sinp(ir — 2d 2 ) — ct^a^ sinp(:r + 2di — 2d 2 )] 

+0(exp(|r|x)), d 2 < x < d 3 , 

p[—a\(X2—ocm-i sinprr - a' l a 2 ...a m ^i s'mp(x - 2di) - ... - aia2—a' m _ 1 
sinp(x — 2d m _i) — a' 1 a! 2 az...a rn -\ smp(x + 2d\ — 2d 2 ) — ... 



-ai...a' i ...a'j...a m -i sinp(a; + 2di — 2dj) 
-ai...a' i ...a' j ...a' k ...a m - 1 sinp(x — 2di + 2dj — 2d k ) + ... 
-oi x oi 2 ...oi m _ x sinp(.x + 2(-l) m - 1 d 1 + 2{-l) m - 2 d 2 + .. 



2d m -t)} 



-0(exp(\T\x)), 



d m - X < X < IT, 



and a, = 



where 

(3.3) oii = 
\ ) 2 '2 

for i — 1,2,..., to— 1. The characteristic function satisfies 

A(A) ~pw(Tr)[aia 2 ...a m ^i sinp7r + a' 1 a 2 ...a m -i sinp(7r — 2di) - 
sinp(7r — 2d m -i) + a' 1 a' 2 a^...a m -i sinp(7r + 2d\ — 2d 2 ) 
+ ai...a' i ...a'j...a m -i sinp(7r + 2d. t - 2dj) 
+ ai...a' i ...a'j...a' k ...a m -i sinp(7r — 2di + 2dj — 2d k ) + ■■ 
+ ai4...a^_! sinp(7r + 2(-l) m - 1 d 1 + 2(-l) m - 2 d 2 + . 

(3.4) +0(exp(|r|7r)). 



a\a 2 ...a r 



2d„ 



Proof. Suppose C(x, A) and S(x, A) are the cosine and sine-type solutions of (|2.1 
with jump conditions (12.31) corresponding to the initial conditions 



(7(0, A) = 1, C"(0, A) = and S(0, A) = 0, S"(0, A) = 1. 

First of all observe 

' cos px + 0( c -^^), 0<x<di, 
a\C\(di,\) cosp(x — di) + -^C[(di,X) sinp(x — d\) 
+0( exp|T| f- dl) ), d 1<x <d 2 , 



a 2 C 2 (d 2 ,X) cos p(x - d 2 ) + ^C' 2 (d 2l \)su\p(x - d 2 ) 



C{x,X) = < 



-0( 



cxp |t|(z— d 2 ) 



d 2 < x < d^ 



a m _iC m _i(d m _i, A) cos p(x - d m -i)+ 



+- s J ±c 'm-i(dm-i, A) sin p(x - d m -\) 

I Q/ gxp]r|(g-dm-i) \ 



d 



m — 1 



< X < 7T. 
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Next we substitute the ?'th statement into the (i + l)'th statement to obtain 

'cospx + 0{^j^), 0<x<di, 

a\ cosp(x) + a[ cosp(x — 2d\) + 0( exp l T l x ) ; d\ < x < d 2 . 
ol\oli cos p(x) + a\a 2 cos p(x — 1d\) + aia 2 cos p(x — 2d 2 ) 
+a[a 2 cosp(x + 2d\ — 2d 2 ) + O ( cxp|T|a " ) , d 2 < x < d 3 , 



C(x,X) = < 



a.\a 2 ...u m -\ cos px + ai...a^...a m -i cosp(a; — 2di) 

2dn 



2d ri 



l) 



+ai...a' i ...a' j ...a m -i cosp(x + 2di - *u, 3J 
+ai...a' i ...a'j...a' k ...a m -i cosp(x — 2dj + 2dj — 2d k ) ■ 
+a' 1 a' 2 ...a' m _ 1 cos p{x + 2{-l) m - 1 d 1 + 2{-l) m ~ 2 d 2 
+0( £S ^), d m -i<x<^, 

where a* and a[ is defined in (|3.3[) and i < j < k, i,j, k — 1,2, ..., m — 1. Similar 
calculations establish the asymptotic form of C'(x, A), S(a;, A), and S'(x, A). This 
proves the theorem upon observing y(x, A) = C(x, A) + /i S(x, A). D 

It follows from the above theorem that 

(3.5) \cp^(x,X)\ =0(|p| l/ exp(|r|x)), < x < tt, f = 0, 1. 

By changing a; to 7r — x one can obtain the asymptotic form of ip(x, A) and ij)'(x, A). 
In particular, 

(3.6) \^ v \x,X)\ =0(\p\ v &cp(\r\(ir-x))), < a: < tt, i/ = 0,1. 
As a consequence of Valiron's theorem ([101 Thm. 13.4]) we obtain: 

Theorem 3.2. The eigenvalues X n = p\ of the boundary value problem L satisfy 

p n = n + o(n) 

as n — > oo. 

4. Uniqueness results for Robin boundary conditions 

In this section we investigate the inverse problem of the reconstruction of a 
boundary value problem L from its spectral characteristics. We consider three 
statements of the inverse problem of the reconstruction of the boundary- value prob- 
lem L : from the Weyl function, from the spectral data {A n ,7„} n >o, and from two 
spectra {X n ,p, n } n >o- 

The Weyl m-function is defined by 



(4.1) 



m(X) 



A(A) 



By (|2.8[) and (|3.6[) we obtain the asymptotic expansion 

1 



(4.2) 



m(A) 



^A 



+ o(\- 1 ) 



along any ray except the positive real axis. 

Let x( x > A) be a solution of ()2.1[) subject to the initial conditions 

X (0,A) = 0, X '(0,A) = 1 
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and the jump conditions (|2.3[) . It is clear that W(cp,x) =1^0 and the function 
ip(x, A) can be represented as 

(4.3) 6(x, A) := ^^ = x (x, A) - m{X) v {x, A). 

The functions 9{x, A) and m(A) are called the Weyl solution and the Weyl function, 
respectively for the boundary value problem L. Clearly 

(4.4) W(cp(x,X),e(x,X))=l. 

Lemma 4.1. The Weyl function m(A) is a meromorphic Herglotz-Nevanlinna 
function, 

(4.5) Im(m(A))=Im(A)||0(A)||^ 
and can be represented as 

oo 

(4-6) m(A)^-^-, 

n— 

where 

(4-7) V - ^ , < oo, V 7 > i 

Proof. The first relation follows after a straightforward calculation using 



(4.8) Im(6»(Tr,A)6l'(7r,A))-Im(6'(0,A)6''(0,A)) = Im(A) / |0(x, A)| 2 w(a;)dx. 

Jo 

Hence m(z) is a Herglotz-Nevanlinna function (i.e. it maps the upper half plane to 
the upper half plane) and by the asymptotics (|4.2j) it has a representation of the 
form (VJ„ Lem. 9.20]) 



m(A) = 
where p is a Borel measure satisfying 



dp(t) 

X n — t 



dp{t) v 7 > i. 



1 + |A|t' ' 2' 

Since by (14.11) the Weyl function is meromorphic it follows that p is a pure point 
measure supported at the poles with masses given by the negative residues. Hence 
the result follows from Lemma [2T2] Q 



Now we are ready to prove our main uniqueness theorem for the solutions of the 
problems (|2.1j) - (|2.3j) . For this purpose we agree that together with L we consider 
a boundary value problem L of the same form but with different coefficients q{x), 
h, H , hi, bi, Ci, di. If a certain symbol r\ denotes an object related to L, then i) will 
denote the analogous object related to L. 

Theorem 4.2. If m(A) = m(A) and w(x) = w(x) then L = L. Thus, the specifi- 
cation of the Weyl function and the weight function w(x) uniquely determines the 
operator. 
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Proof. It follows from flU) and ((43J that 

(4.9) |^)( 2 ;,A)|<C|pr- 1 exp(-|r|a ; ), 1/ = 0, 1, 

as A — ► oo along any ray except the positive real axis. Define the matrix P(x, A) 
[Pjk(x, \)]j,k=i,2 by the formula 

P(rr \\ ^( X ' A ) H*,X)\ _ /V(z,A) 6»(a;,A) 
n ' J ^(a:,A) 0'(z,A)J vVM 0'Oz.A) 

Taking (|4.4[) into account we calculate 

r Pn{x,X) P 12 (x,X)\ _ fipd'-ifi'd <p6-<p§ 



(4 ' 10) V p 2i(^,A) P 2 2{x,X)J \<p'6'-Cp'6' <p6'-ip' 

and 

'(fix, A)\ _ (Pn(x, \)<p(x, A) +Pi 2 (x, \)<p'(x, A) 



(4-i1 -' l ^(^,A); V- p ii( a: > A ) fl ( a: > A ) + p i2(a ; .^(a : ^) 

It is easy to see that the functions Pjk{x, A), j, k = 1, 2, are meromorphic in A with 
simple poles in the points A ra and A n . Moreover, if 771(A) = fh(X), then from (|4.3|) 
and (I4.10J) . Pn(a:, A) and Pi2(£, A) are entire functions of growth order 1/2 in A. 
From (HU) 

(4.12) |Pii(!t,A)|<C, |Pi 2 (x,A)|<^ 

along any ray except the positive real axis. Moreover, by our hypothesis this func- 
tion has an order of growth s and thus we can apply the Phragmen-Lindelof theo- 
rem (e.g., |10| Sect. 6.1]) the two half-planes bounded by the imaginary axis. This 
shows that the functions Pn and P12 are bounded on all of C and thus constant by 
Liouville's theorem. Since P12 vanishes along a ray it must be zero and we obtain 
Pu(x, A) = A(x) and Pi 2 (x, A) = 0. Using (|4TT|) . we get 

(4.13) ip(x, A) = A(x)tp(x, A), 9{x, A) = A(x)0(x, A). 

It follows from (12TTT1) . W(tp(x,X),6(x,X)) = W(tp(x,X),9(x,X)) = 1 and so we 
deduce A(x) = "^X = 1, that is, tp(x. A) = (p(x, A), 9(x, A) = 9(x, A), and ip(x, A) = 
■4>(x,X). Therefore from (|2TJ) . (|2~3]) . (l2~TTj) , and ([121) we get g(a;) = ?(»), a.e. on 
[0, 7r] and a, = a^, &i = 6j, Cj = Cj, and g^ = di for i = 1,2, ...,m — 1, /i = h and 
H = H for j = 1, 2, 3. Consequently L = L. □ 

Note that this theorem is optimal in the sense that the weight function cannot be 
determined from 771(A) since a unitary transformation as in Lemma 12.31 can be used 
to change the weight without changing tti(A). Note that the condition w(x) = w(x) 
will hold if we have for example a^ = dibi = 1 for all i or aibi = dibi, and di = di, 
for I = 1,2, ...,77i— 1. 

By virtue of Lemma T4. II we also get: 

Corollary 4.3. If X n = X n and j n = j n , for n = 0, 1, 2, ..., and w(x) = w(x) then 
L = L. 
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Finally, let us consider the boundary value problem L k which is the problem 
where the boundary condition L\(y) is replaced by 

vi )= jV(°) + fc y(o) = o, fceK > 

lW> \»(0) = 0, fc = oo. 

Let {/J, n }n>o be the eigenvalues of the problem L k . 

Corollary 4.4. Suppose k ^ h. If A n = A„ and fi n — fi n for n = 0,1,2,..., and 
w(x) — w(x), then L = L. 

Proof. We begin with the case k = oo. The numbers A„, fi n are the poles and 
zeros of m(A) and hence determine it uniquely up to a constant by Krein's theorem 
[TU1 Thm. 27.2.1]. This unknown constant can be determined from f|4.2[) . The case 
k ^ h follows in the same manner using m(A) + (fc — h) . □ 

Finally, we are also able to extend Hald's theorem to the case of finitely many 
transmission conditions. 

Theorem 4.5. If \ n = X n , w(x) = w(x), L\ = L\, q(x) — q(x) for a.e. x < ^ 
and Ui = Ui, Vi = Vi for all i with di < ^, then L = L. 

Proof. By the Hadamard factorization theorem W((p, ip) = C W(<p, tp) for some 
constant C which can be determined from the asymptotic as A — > oo: 



c= TT ^>o. 



Oii 
i:di>7r/2 



Furthermore, our assumptions imply 



71 



for some entire function F(X) of growth order at most i. Solving for F and taking 
the limit x f § we obtain 

^(|-,A)-CV>(f-,A) VX|-,A) / ^(f- A) _ \ 

1 j </?(§- A) " ^(|-,A) \C^-,X) )■ 

Now the expression in parenthesis vanishes along every ray different from the pos- 
itive real axis while the expression in front remains bounded by the asymptotics 
(|3.ip for tp and the analogous result for ip, ip. Thus it must be identically zero by the 
Phragmen-Lindelof theorem. Finally, x(x, A) = x(x, A) for x < ? implies that the 
associated Weyl functions are equal and the claim follows from Theorem 14.21 □ 

5. Uniqueness results for eigenparameter dependent boundary 

conditions 

In this last section we will replace the Robin boundary condition (12.21) by the 
following eigenparameter dependent boundary conditions 

Li(i/) := A(y'(0) + h iy (0)) - h 2 y'(0) - h 3 y(0) = 0, 

(5.1) L 2 (y) := X(y'(7f) + H iy (n)) - H 2 y'{-K) - H 3 y(n) = 0, 
where we assume that hj, Hj, j = 1, 2, 3, are real numbers, satisfying 

(5.2) ri := h 3 - h x h 2 > and r 2 := H X H 2 - H 3 > 0. 



10 



M. SHAHRIARI, A. JODAYREE AKBARFAM. AND G. TESCHL 



In order to obtain a self-adjoint problem we will use the following Hilbert space 
H := ^2((0, 7r); w) © C 2 with inner product defined by 

(5.3) 

<f{x)\ (g(x? 

r = I A , G = 51 

h) V -92 



(^G) 



H 



/fi™ 



w(0) _ w(tt) _ 
-Ji9i H J2ff2, 



n 



r-2 



1/2 

Again the associated norm will be denoted by \\F\\-u — (F^ F)^ . Next we introduce 

Ri(y) ■= i/(o) + fci»(o), ^i(y) := W(o) + ^y(o), 

fl 2 (j/) := i/'(tt) + JJu/fTr), #,(</) := F 22 /(?0 + ffal/M- 
In this Hilbert space we construct the operator 



(5.4) 

with domain 
(5.5) 



dom (A) = { F = 




by 



if 
AF = | R[(f) 

,W)j 



A-.n^n 



f,f € ^(Uo" 1 - 1 (dj.di+i)), £/ € £ 2 (0,vr) 

£/*(/) = ^(/) = o, A = Ri(f), h = Mf) 



fix) 

with F = | i?i(/) | e dom (A) . 

Mf)j 



By construction, the eigenvalue problem for A, 
(5.6) AF = AF, Y :- 



y{x) 

R^y) | Gdom(A), 
^a(l/)> 



is equivalent to the eigenvalue problem (|2.1j) . (|2.3[) . and (|5.1[) for L. A straightfor- 
ward calculation shows: 

Lemma 5.1. The operator A is symmetric. 

In particular, the eigenvalues of A, and hence of L, are real and simple. To 
see that they are simple it suffices to observe that the associated Cauchy problem 
(|2~Tj) . ([2~3]) subject to the initial conditions f(x ± 0) = f , f(x ± 0) = f x (with 
.To € [0, 1]) has a unique solution. 

Suppose that the functions ip(x, A) and i/j(x, A) are solutions of (|2.1j) under the 
initial conditions 

(5.7) p(0, A) = A - h 2 , cp'(0, X) = h 3 - Xh u 
and 

(5.8) ^(tt, \)=H 2 - A, V'(tt, A) = Ai/i - H 3 

as well as the jump conditions (|2.3p . respectively. Moreover, we set 

(5.9) A(A) := W(¥>(A),V(A)) = »ijW)). 
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Then A(A) is an entire function whose roots A„ coincide with the eigenvalues of 
L. Moreover, the eigenfunctions ifi(x, X n ) and ipi(x,X n ) associated with a certain 
eigenvalue A„, satisfy the relation ipi{x,X n ) = /3 n ipi(x,X n ), where, by (|5.7[) . 



(5.10) 



V>'(0,A n ) + M>(0,A n ) 



n 



We also define the norming constant by 



7 „:=||$(x,A n )|| 



$(z,A)= ( R^ip) 



Then it is straightforward to verify: 



Lemma 5.2. All zeros X n of A(A) are simple and the derivative is given by 



(5.11) 



A(a„) = - 7l ;X- 



The same argument as for Theorem 13. II shows: 



Theorem 5.3. Let A — p and t := Imp. For equation (|2.1[) with boundary 
conditions (|5.ip and jump conditions (|2.3|) as |A| — > 00, t/ie following asymptotic 
formulas hold: 



(5.12) 



p 2 cospx + 0(pexp(\r\x)), 



< x < di , 



p 2 [ai cospx + «j cosp(x — 2di)\ + 0(pexp(|r|a;)), 



di < a; < <i2 



p 2 [aia2 cospx + 0^0:2 cosp(x — 2di) + a\a' 2 cos p(x — 2&2) 



<p(x,A) = < 



-o! x q!<2 cos p(x + 2<ii — 2rf 2 )] + 0(pexp(|r|x)), 



p 2 [aia2...a m _i cospx+ 

+a' 1 a2---a»n-i cosp(x — 2di) + ... 
+a 1 ct2...a' m _ 1 cosp(x — 2d m _i)+ 
+a' 1 a 2 a3...a m _i cosp(x + 2di — 2g?2) 



d 2 < x < d 3 , 



-a\...a. 



.a m _i cosp(x + 2di — 2dj 



-ai...a' l ...a' J ...a' k ...a m -i cosp(x — 2d % + 2dj - 2d k ) + ... 
-aia 2 ...a^_! cosp(x + 2(-l)'"- 1 d 1 + 2(-l) m - 2 d 2 + ... - 2d m _i)] 



+0(pexp(|r|x)), 



d m _i < x < IT, 
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and 
(5.13) 



ip'(x,X) 



< x < di , 

o2, 



p 3 [— sin px] + 0(p 2 exp(|r|x)), 

p 3 [— ai sin pa; — a[ sinp(a; — 2di)] + 0(p 2 exp(|r|a;)), d\ < x < d 2 , 

p 3 [— ot\OL2 sin pa; — a[a 2 sinp(a; — 2d\) — 

—a\a' 2 sinp(a; — 2d 2 ) — ct[a 2 sinp(a; + 2d\ — 2d 2 )] 

+0(p 2 exp(|r|a;)), d 2 < x < d 3 , 



p 3 [—aia 2 ...a m ~i sin pa; — a' 1 a 2 ...a m -i sinp(a; — 2d\) — ... — aia 2 ...a' m 
sinp(a; — 2<i m _i) — a l 1 cJ 2 az...a. rn -\ sinp(a; + 2d\ — 2d 2 ) — ... 
— ai...a' i ...a'j...a m -i sinp(a; + 2di — 2dj) 
—a>i...a' i ...a'j...a / k ...a m -i sinp(a; — 2di + 2dj — 2dk) + ... 
-a x a 2 ...a' m _ x sinp(a; + 2(-l) m - 1 d 1 + 2{-l) m - 2 d 2 + ... - 2d m -x)] 



+0(p 2 exp(|r|x)), d m _i < x < w, 



where 
(5.14) 



a-i + b t 



and a'; 



2 , 2 , 

for i — 1, 2, ..., m — 1. The characteristic function satisfies 
A(A) =p 5 w(ir)[aia 2 ...a m -i sinp7r + a 1 a 2 ...a m -i sinp(7r — 2d\) 
sinp(7r — 2d r „_i) + a l 1 a' 2 a 3 ...a m -i sinp(7r + 2d x — 2d 2 ) - 
+ a\...ct i ...cty..a m -\ sinp(7r + 2di — 2dj) 
+ ai...a' i ...a'j...a' k ...a m -i sinp(7r — 2d{ + 2dj — 2dk) + ... 



a±a 2 ...a 



m — l 



«1"2- 



m—l 



sinp(7r + 2(-l) m - i d 1 + 2(-l)" l ^ 2 + ... - 2d m _i) 



(5.15) 



0(p 4 exp(|r|7r)). 



It follows from the above theorem that 

(5.16) |<p M (a;,A)| =0(|pr +2 exp(|r|a;)), < x < n, v = 0, 1 

and so by substituting x with n — x we get the asymptotic form of ip{x, X) and 
^'(x,X). In particular, 

(5.17) |V> (v) (a;,A)| = 0(|pr +2 exp(|r|(7r - a;))), < x < n, v = 0, 1. 
As a consequence of Valiron's theorem ([TQl Thm. 13.4]) we obtain: 

Theorem 5.4. The eigenvalues X n — p 2 of the boundary value problem L satisfy 

p„ = n + o(n) 
as n — > 00. 

The Weyl rn-function is defined by 

Ri(tp(X)) ^'(0,A) + M>(0,A) 



(5.18) 



m(A) = -- 



riA(A) 



riA(A) 
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By (12. 8[) and (|5.17|) we obtain the asymptotic expansion 

(5.19) m(A) = — L + 0(A^ 3 / 2 ) 

r\\ 

along any ray except the positive real axis. 

Let x( x i ^) be a solution of (|2.ip subject to the initial conditions 

x(0,A) = -i, x '(0,A) = ^i 
r\ n 

and the jump conditions (|2.3p . It is clear that W(ip,x) = 1^0 and the function 
i>(x, A) can be represented as 

(5.20) 0(x, A) := ^^ = x (z, A) - m(\)<p(x, A). 

The functions #(x, A) and m(X) are called the Weyl solution and the Weyl function, 
respectively for the boundary value problem L. Clearly 

(5.21) W(cp(x,\),6(x,X)) = l. 

Lemma 5.5. The Weyl function m(A) is a meromorphic Herglotz-Nevanlinna 
function, 

I e(x,x) \ 

(5.22) Im(m(A)) = Im(A)||e(A)||^, @(x,\) = Ri(6(X)) , 



MOW)/ 



and can be represented as 



(5.23) m(A) = £_*!_ 



71=0 



where 

oo 1 

(5.24) E7n = ~. 

n=0 

Proof. The first relation follows after a straightforward calculation using 



(5.25) Im(0(7r,A)0'(7r,A))-Im(0(O,A)0'(O,A)) = Im(A) / \6(x, \)\ 2 w{x)dx. 

Jo 

The first part follows as in the proof of Lemma 14.11 Computing the asymptotic of 
(T5T23T) and comparing with (f5"7ft)|) shows ([ST24)) . □ 

Now we are ready to prove our main uniqueness theorem for the solutions of 
the problems (12. ip . (|2.3[) . and (|5.1[) . For this purpose we agree that together with 
L we consider a boundary value problem L of the same form but with different 
coefficients q(x), hj, Hj, hi, bi, Ci, di. 

Theorem 5.6. // m(A) = ifi(X) and w(x) — w(x) then L = L. Thus, the specifi- 
cation of the Weyl function and the weight function w{x) uniquely determines the 
operator. 

Corollary 5.7. // A ?l = A n and j n = j n , for n = 0, 1, 2, ..., and w(x) = w(x) then 
L = L. 
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Finally, let us consider the boundary value problem L' , where we take the con- 
dition y'(0) + hiy(0) = instead of the condition (|5.1j) in L. Let {/x„} n >o be the 
eigenvalues of the problem V . 

Corollary 5.8. If X n — X n and fi n = \x n for n = 0,1,2,..., and w(x) = w(x), 
7'i = f\ then L = L. 

Theorem 5.9. If X n = A„, w(x) = w(x), L\ = L\, q(x) = q(x) for a.e. x < ^ 
and Ut = Ui, Vi = Vi for all i with d t < ^, then L = L. 
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Addendum 



In the introduction we attribute Theorem l4.2l and its two corollaries in the special 
case for one transmission condition with determinant one (i.e., w = 1) to [3]. After 
publication of this paper we learned that this special case was already obtained 
earlier by Yurko in 

• V. A. Yurko, Integral transforms connected with discontinuous boundary 
value problems, Integral Transforms Spec. Funct. 10, 141-164 (2000). 
We regret we only learned about this paper after publication of our article. In 
addition, we note that Yurko's paper also solves the inverse spectral problem in 
this context. 
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